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Abstract 

The dipole-coupled two-level atoms(qubits) in a single-mode resonant cavity is studied by ex- 
tended bosonic coherent states. The numerically exact solution is presented. For finite systems, 
the first-order quantum phase transitions occur at the strong interatomic interaction. Similar to 
the original Dicke model, this system exhibits a second-order quantum phase transition from the 
normal to the superradiant phases. Finite-size scaling for several observables, such as the average 
fidelity susceptibility, the order parameter, and concurrence are performed for different interatomic 
interactions. The obtained scaling exponents suggest that interatomic interactions do not change 
the universality class. 

PACS numbers: 42.50.Nn, 64.70.Tg, 03.65.Ud 
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I. INTRODUCTION 



The coherent emission from under damped one- dimensional Josephson arrays coupled to 
a single-mode electromagnetic cavity has been experimentally studied[l|. This system re- 
sembles that of the Dicke model |2j, which describes N two- level atoms (qubits) coupled to 
a cavity field. It has been shown |3J that a modified Dicke Hamiltonian including a dipole- 
dipole interaction between the junctions can describe the cavity-junction system better than 
the original Dicke model. 

The modified Dicke Hamiltonian is also an extension of effective two-qubit modelQ]. A 
small number of coupled qubits in a single-mode cavity is also very interesting, because it 
may be realized in several solid-state systems, such as an ensemble of quantum dots Bose- 
Einstein condensates jsj, coupled arrays of optical cavities used to simulate and study the 



behavior of strongly correlated systems [Tj], and the superconducting quantum interference 
device coupled with a nanomechanical resonator g), 0] . 

The original Dicke model without a rotating-wave approximation (RWA) on the large 



scale was exactly solved by the present authors in the numerically sense [10|. Most works 
on the modified Dicke model are limited to the rotating-wave approximation. With the 
progress of the fabrication, the artificial atoms may interact very strongly with on-chip 
resonant circuits . Il2|. the rotating-wave approximation can not describe well the 

strong coupling regime 13], so the numerically exact solution to the modified Dicke model 



without rotating-wave approximation is also of considerable significance and highly called 
for. 

Quantum phase transitions (QPTs) in the original Dicke model has attracted considerable 



attentions recently 



1J-|20|. With the consideration of additional interatomic interactions, one 



nature question is its effect on the modified Dicke model. 

In this paper, we extend our previous exact technique to solve the finite size modified 
Dicke model where dipole-dipole interaction between the qubits are taken into account. The 
QPTs are then studied systematically. The paper is organized as follows. In Sec. II, the 
numerically exact solution to the finite-size modified Dicke model is proposed in detail, and 
the analytical solution in the thermodynamic limit is also presented. The numerical results 
for both small and large system size are given in Sec. Ill, where the characterization of the 
QPTs also also performed. The brief summary is presented finally in the last section. 
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II. MODEL HAMILTON! AN 

dipole-coupled N two- level atoms interacting with a single-mode cavity Q, Q] can be de- 
scribed by the following modified Dicke model [3] 

A N A N Q N 

H = uja+a + - J] a\ + -= (a+ + a) £ o» m + — + , W 



2 f-f 2 v^V v ; f-f x 27V . 

i=i v i=i 

where a + and a are the field annihilation and creation operators, A and u are the transition 
frequency of the qubit and the frequency of the single bosonic mode, A is the coupling 
constant of the atom and cavity, Q is the interacting strength of two two-level atoms, and 
o % k {k = x,y, z) is the Pauli matrix of the i junction. For convenience, the Hamiltonian can 
be rewritten in terms of the collective spin operators: Sk = 5^ i=1 cr|./2, (k = x,y,z) and 
S± = S x i iSy 

2 A 90 
H = ua + a + AS Z + -= (a+ + a) S x + — (S 2 - S 2 Z ), (2) 
y N A 



Where a constant is neglected. Without atom-cavity coupling, i.e. A = 0, we have 

2Q 2 2 

Hatoms = AS Z + — (S — S z ), (3) 

Note that the Dicke state {\j, m) ,m — —j, — j + 1, ...j — 1, j} is the eigenstate of S 2 and S z 
with the eigenvalues j(j + 1) and m. The eigen energy of H atoms is given by 

E = Am + —(3(3 + 1) -m 2 ), (4) 

The ground-state energy is easily obtained as 

So in the ground-state, we have j = for ^- > A and j = iV/2 otherwise. Actually, Eq. (3) 



is just the Hamiltonian of the isotropic Lipkin-Meshkov-Glick model 2l|, where QPT is of 
the first-order, as is clearly shown above. 



Next, we use a transformed Hamiltonian with a rotation around an y axis by an angle 
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, so S x — y S Z} S z — y — S x the Hamiltonian now reads 



H> = coa+a - ^(S + + 8.) + (a+ + a) S z + ^(S 2 - S 2 ), (5) 



where S + and £_ are the collective spin raising and lowing operators and obey the SU(2) Lie 
algebra [S + , £_] = 25*2, [5*2, S±] = ±S±. Since [H 1 , S 2 ] = 0, we suppose S 2 = j(j + 1), where 



the total angular momentum j = N/2 — r. Due to the interaction between any two-level 
atoms, similar to the above discussion for A = 0, we have no reason to set r = 0, unlike in 



the Dicke model[l0|. In the practical calculations for the ground-state, r is regarded as a 
variational integer number and will be determined by the minimization of the ground-state 
energy. So the Hilbert space of this algebra is spanned by the Dicke state { \j, m) , to = 
—j, —j + I, ...j — 1, j}, which is the eigenstate of S 2 and S z with the eigenvalues j(j + 1) 
and rn. 

The Hilbert space of the total system can be expressed in terms of the basis 
{|^m)b *S> \h m )}i where the state of bosons \(f m ) b and the integer j are to be determined. A 
"natural" basis for bosons is Fock state \l) b = [{a + ) l /\/l\} |0) 6 . As in the Dicke model, the 
bosonic number here is not conserved either, so the bosonic Fock space has infinite dimen- 
sions, the standard diagonalization procedure (see, for example, Ref. [ljj] ) is to apply a 
truncation procedure considering only a truncated number of bosons. Typically, the conver- 
gence is assumed to be achieved if the ground-state energy is determined within a very small 
relative errors. Within this method, one has to diagonalize very large, sparse Hamiltonian 
in strong coupling regime and/or in adiabatic regime. Furthermore, the calculation becomes 
prohibitive for larger system size since the convergence of the ground-state energy is very 
slow. Interestingly, this problem can be circumvented in the following procedure. 

By the displacement transformation A m = a+g m with g m = 2\m/u\/N~, the Schro dinger 
equation can be described in columnar matrix, and its m row reads 

\(f m ) b \j, m - 1} - Aj+ \(p m ) b \j, to + 1) 
2VL _ _ , . . . 2fi , . , . . 

-J^jmJm-l Wm)b |j» TTl - 2) - —jZ]rn+l Wm) b \j, TO + 2) 

2 

+U (A+A m - g 2 m ) \(f m ) b \j, m) + -N + 1) - {jt-lJm + Jm+lJm)] \<Pm)b b'» m ) 

= E\ip m ) b \j,m) , (6) 



where = |a/(j(j + 1) — m(m ± 1). Left multiplying {(n, j\} gives a set of equations 

2Q 20 

-Aj n \<Pn-l) b ~ Aj+ \¥n+l) b ~ JjrjnJn-l \fn-2) b ~ Jj-jnJn+1 \Vn+2) b 

+u (A+A n - gl) \ip n ) h + -N + 1) - (j+_dn + Jn+iJn )] M b = E \<p n ) b , (7) 

where n = — j, — j + 1, ...j. 

Note that the linear term for the bosonic operator a(a + ) is removed, and a new free 
bosonic field with operator A(A + ) appears. In the next step, we naturally choose the basis 



in terms of this new operator, instead of a(a + ), by which the bosonic state can be expanded 

as 

N tr 



k=0 

Ntr -. 

= E c ^^n( a+ + 9n) k e^ a+ -^ 2 \0) a , (8) 

where N tr is the truncated bosonic number in the Fock space of A(A + ). As we know that the 
vacuum state |0)^ is just a bosonic coherent-state in a(a + ) with an eigenvalue <7 n [10|. So this 
new basis is overcomplete, and actually does not involve any truncation in the Fock space of 
a(a + ), which highlights the present approach. It is also clear that many-body correlations 
for bosons are essentially included in extended coherent states (5). Left multiplying state 
A n (l\ yields 

jw(Z - g 2 n ) + [](] + 1) - {jt-lin + Jn+dn )] J C n ,l 

Ntr 

~ A J2(jn A n (l Ik)^ Cn-l,k + 3n A n (l \k) An+1 C n+1)k J 
k=0 

~ AT S yn3n-\ A n (I \k) An „ 2 Cn-2,k + jt Jn+l An ( l \ k ) A n+2 C n+2,k) = Ec n>h (9) 



fc=0 



where 



with 



A n (l\k) An _ x = (-l) l D Lk (G), An (l\k) An+i = (-l) k D Lk (G), 

An(l\k) An _ 2 = (-l) l D Lk (2G), An (l\k) An+2 = (-l) k D Lk (2G), (10) 



min[lM ' -i)- r Vm.G l+k - 2r „ 2A 



l ' k (1 - rV(b--AM ' 



_ n (I - r)\(k - r)\r\ uj^/N 
Eq. (6) is just a eigenvalue problem, which can be solved by the exact Lanczos diagonal- 
ization approach in dimensions (2j + l)(N tr + l)[ll3] • Note that the eigenvalue problem in 
the pure Dicke model can be reproduced if set Q = 0. As before, to obtain the true exact 
results, in principle, the truncated number N tr should be taken to infinity. Fortunately, it 
is not necessary. It is found that finite terms in state (5) are sufficient to give very accurate 
results with a relative errors less than 10 -6 in the whole parameter space. We believe that 
we have exactly solved this model numerically. 



Thermodynamic limit- In order to study the quantum phase transition explicitly of 
this model, we should evaluate the transition point in the thermodynamical limit iV — > 

00. First, we use the Holstein-Primakoff transition of the collective angular momentum 
operators defined as S + = cfiyN — aJa, S- = y/N — a^aa, and S z = a^a — N/2, where 
[a, a*] = 1. Second, we introduce shifting boson operators and with properly scaled 
auxiliary parameters a and /3 such that = + y/N a and d) = a) — y/~N f3 to describe the 
collective behavior of the Hamiltonian in Eq.( 1). Finally, by means of the boson expansion 
approach, we expand the H with respect to the new operators and d) as power series in 
1/N. According to Hamiltonian (2), we have the scaled ground state energy 

^1 = ua 2 - AXaPy/l^W + A(/3 2 - i) - 2fi(/3 2 - \f + 5 

The critical points can be determined from the equilibrium condition d[Eo(a, (3)/N]/da = 
and d[E (a, (3)/N]/d(3 = 0, which leads to two equations 

ua - 2A/V1 - I3 2 = 0, 

2aA v / l -P 2 - 2a\(3 2 - f3A + 4Q(3(f3 2 - -) = 0. (11) 

\Jl-(3 2 2 

Then we can obtain the critical atom-cavity coupling constant in the second-order QPT 
A c = y/u(A + 2Q)/2 

III. RESULTS AND DISCUSSIONS 

The exactly numerical results are presented in this section to study the properties of QPT 
in this model. Without loss of generality, we set u = A in the whole calculation. The units 
are taken of u = 1 for convenience. 

First, we study a finite number of coupled two- level atoms in a cavity. The maximum 
value of j is j max = A/2. It is very surprising that j ^ j max in the ground-state for finite 
systems, unlike in the pure Dicke model. In Fig. [TJJa) and (b), we plot the value of j as 
a function of A and Q' = 2Q/N in the ground-state for A = 4. In the weak interaction 
of the two- level atoms with small Q', j = jmax-, similar to the pure dicke model where the 
interaction of the two- level atoms is neglected. As Q' increases, the value of j is reduced. In 
the intermediate coupling range of atom-cavity A, the value of j is reduced gradually by step 

1, until to zero in the strong atomic interaction. This phenomena disappears in the strong 
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FIG. 1: The total momentum angular in the ground-state j as a function of the coupling constant 
A and the interaction strength f2 between atoms for N = 4. 
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FIG. 2: The average ground-state energy E/N, its second derivatives d 2 E/d\ 2 , and the total 



angular momentum j as a function of A for different atomic interactions Q = 2.2 (a), f2 = 2.5 (b), 
and £1 = 3.0 (c) for N = 4. The dashed lines denote the energy for the total angular momentum 
j = jmax = N/2, which is clearly higher the true ground-state energy. 



coupling range of atom-cavity A, where j is always equal to j max - For A = 0, according to 
Eq. (4), the value of j jumps to zero when the 2Q/N = 1, consistent with the observation 
in Fig. [TJ The above observations hold true for any finite number of atoms. 

The jump of j in the A — Q plane exhibited in Fig. [T] may be related to the first-order 
QPT in a finite system. It is known that there should be a second-order QPT in this systems 
without interatomic coupling. Since both the first- and second-order QPTs can be simply 
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characterized by singularities of the ground-state energy, we calculate ground-state energy 
and its second derivative. These results together with the values of j are collected in Fig. 
[5] for three typical value of Q above the first-order QPT point at A = 0. For smaller vale 
of Q = 2.2, two jumps of j indicates two first-order QPTs, as shown in Fig. Ufa). The 
non-analytic feature of the ground-state energy and sudden drop of its second derivative at 
the two jump points give the evidence of the first-order QPT. Above the last jump point, the 
ground-state energy are continuous and its second derivative shows a smooth drop around 
A c , demonstrating a sign of the second-order QPT. \ C (Q, N) is also presented in Fig. [TJ 
which divides the j = j max regime into two parts. However, as shown in Fig. [2]^c), for 
larger value of Q = 3, there is only one jump of j from to j ma x — N/2 in the whole 
coupling regime. Both non-analytic feature of the ground-state energy and sudden drop of 
its second derivative suggest a first-order QPT at this jump point. No sign of second-order 
QPT is observed in this case. For the intermediate vale Q = 2.5, two first-order QPTs 
occur sequentially with the jump of j. We think that both the first and the second-order 
QPT occur in the same point for Q = 2.5 and 3.0, and finally only the characteristic of the 
stronger first-order QPT shows up at the last jump. 

If we naively set the total angular momentum to be j max = N/2 as in the pure Dicke 



model 



101 ]. we can also calculate "ground-state", which are also given in Fig. [2] with dashed 



lines. It is really higher than the true the ground-state energy. The continuous behavior 
in the whole coupling regime is observed in these dashed curves, demonstrating the above 
interesting feature in a finite system stems from the variation of the total angular momentum. 

Next, we will study the effect of the interatomic interaction on the second-order QPT in 
the thermodynamical limit. We focuss on the question whether the interatomic interaction 
alters the universality class of the second-order QPT. 

Recently, the fidelity, a concept in quantum information theory, has been extensively used 
to identify the QPTs in various many-body systems from the perspective of the ground- 
state wave functions 22M29j] . In a mathematical sense, the fidelity is the overlap between 
two ground states where the transition parameters deviate slightly. However, the fidelity 
depends on a arbitrary small amount of the transition parameters, which in turn yields an 
artificial factor. Zanardi et al 1241 introduced the Riemannian metric tensor and You et al 



251 ] proposed the fidelity susceptibility (FS) to avoid this problem independently. 



To analyze the QPT, we first illustrate the scaling behavior of the average FS. The 



S 




FIG. 3: Finite-size scaling of the average FS according to Eq. (fT2|) at the critical point for (a) O = 0, 
(b) n = 0.25, and (c) ft = 0.5. 



finite-size scaling ansatz for the average FS take the form 29] 



^ = f[N"(\ - A max )] (12) 

Xf 

where x™ ax ^ s the va lue of average FS at the maximum point A max , / is the scaling function 
and v is the correlation length critical exponent. This function should be universal for large 
N in the second-order QPTs, which is independent of the order parameter. As shown in Fig. 
[3] an excellent collapse in the critical regime is achieved with the use of v = 2/3 according to 
Eq. (ll2p in the curve for different large size for three values of Q. It is demonstrated that v is 
a universal constant and does not depended on the parameter Q, suggesting the interatomic 
interaction does not change the universality class. 

We then perform the finite size scaling analysis on the order parameter of the QPT, i.e. 
the expectation value of the photon number per atom in the ground-state (cfta)/N. In the 
thermodynamic limit, this quantity changes from zero to finite value smoothly when crossing 
the critical point. In Fig. HI we present this quantity as a function of N for different values 
of Q in log-log scale. Derivatives of these curves are plotted in the inset. The exponent of 
the order parameter is estimated to be —0.66 ± 0.01 for three values of Q, provided another 
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FIG. 4: Scaling of the order parameter (a^a) /N as a function of iV at the critical point for $7 = 0, 
0.25, and 0.5. The inset shows the slope versus 1/N. 

piece of the evidence that the interatomic interaction does not change the universality class. 

Finally, we calculate the scaled concurrence (entanglement) and perform the correspond- 
ing finite size scaling analysis as in Ref. Q, Q, 19 |. The concurrence quantifies the en- 



tanglement between two atoms in the atomic ensemble after tracing out over bosons. In 
the thermodynamic limit, the scaled concurrence at critical point Coo(A c ) can be easily de- 
termined with the solution of Eq. (11). For comparison, we can calculate the quantity 
Coo (^c) — CV(A C ). In Fig. |5l we present this quantity as a function of N for different values 
of Q in log- log scale. Derivative of these curves is presented in the inset, and the exponent of 
concurrence is estimated to be —0.33 ± 0.005 for three values of Q. This again demonstrates 
that the universality class is not altered by the interatomic interaction. 



IV. CONCLUSION 



In summary, a finite number of dipole-coupled two-level atoms in a single-mode resonant 
cavity is solved exactly with the use of extended bosonic coherent in the numerically sense. 
A number ( the maximum value is N/2) of first-order QPTs occur in this system with N 
coupled atoms, different from in the pure Dicke model. The total angular momentum in the 
ground-state is altered with the coupling constant and the interatomic interactions. The 
quantum criticality is also studied in terms of the ground-state fidelity, order parameter 
and concurrence in very large systems up to 2 10 or more. The finite-size scaling analysis 
for these observables are performed. The corresponding scaling exponents obtained remains 
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FIG. 5: Scaling of concurrence as a function of N at the critical point for 17 = 0, 0.25, and 0.5. 
The inset show the slope versus 1/N. 

unchanged with the interatomic integrations, demonstrating that the university class is not 
changed. It should be pointed out that all eigenfunctions and eigenvalues obtained in the 
modified Dicke model on the large scale might also be used to explore the mechanism for 
the coherent radiation in one-dimensional Josephson arrays coupled to a single-mode cavity 
at both zero and finite temperatures, which may be our future work. 
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